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The Lagrangian, the Hamiltonian and the constant of motion of the gravita- 
, Ph | tional attraction of two bodies when one of them has variable mass is con- 

sidered. This is done by choosing the reference system in one of the bodies 
which allows to reduce the system of equations to 1-D problem. The trajecto- 
ries found in the space position- velocity, (x, v), are qualitatively different from 
^ ! those on the space position-momentum, 
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1. Introduction 



Mass variable systems has been important since the foundation of the classical 
mechanics and have been relevant in modern physics [1]. Among these type 
of systems one could mentioned: the motion of rockets [2], the kinetic theory 
of dusty plasmas [3], propagation of electromagnetic waves in dispersive and 
nonlinear media [4], neutrinos mass oscillations [5], black holes formation [6], 
and comets interacting with solar wind [7]. The interest in this last system 
comes from the concern about to determinate correctly the trajectory of the 
comet as its mass is changing. This system belong to the so called two-bodies 
problem. The gravitational two-bodies system is one of the must well known 
systems in classical mechanics [8] and is the system which made a revolution 
in our planetary and cosmological concepts. Comets loose part of their mass 
as traveling around the sun (or other star) due to their interaction with the 
solar wind which blows off particles from their surfaces. In fact, it is possible 
that the comet could disappear as it approaches to the sun [9]. So, one must 
consider the problem of having one body with variable mass during its gravita- 
tional interaction with other body. This problem has received many attention 
and is known as Gylden-Meshcherskii problem [10]. However, little or non at- 
tention at all has been made about the possible qualitative different picture of 
the trajectories on the space position- velocity and position-momentum spaces 
generated by the constant of motion and Hamiltonian of the system which, we 
think, may have interest in relation with the study of consistent Hamiltonians 
for autonomous systems. 

In this paper, one considers the problem of finding the constant of motion, 
Lagrangian, and Hamiltonian, for the gravitational interaction of two bodies 
when one of them is loosing its mass during the gravitational interaction. The 
mass of one of the bodies is assumed much larger than the mass of the other 
body. Choosing the reference system on big-mass body, the three-dimensional 
two-bodies problem is reduced to a one-dimensional problem. Then, one uses 
the constant of motion approach [11] to find the Lagrangian and the Hamil- 
tonian of the system. A model for the mass variation is given for an explicit 
illustration of these quantities. With this model, one shows that the trajec- 
tories in the space position-velocity (defined by the constant of motion) are 
different than the trajectories on the space position-momentum (defined by 
the Hamiltonian). 
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2. Reference system and constant of motion 

Newton's equations of motion for two bodies interacting gravitationally, seen 
from arbitrary inertial reference system, are given by 

d ( drA Gm 1 m 2 . . . , 

^r^j = -R-^ (ri - r2) (ia) 

and 

d ( dr 2 \ Gm 1 m 2 

^r 2 ^J = -^7F (r2 - ri) ' {lb) 

where mi and m 2 are the masses of the bodies, r x = (xi,yi,Zi) and 
r 2 — (^2,2/2,^2) are the vectors position of the two bodies from our reference 
system, G is the gravitational constant (G = 6.67 x 10~ u m 3 /Kg sec), and 
I ri — r 2 | is the Euclidean distance between the two bodies, 



ri - r 2 



= |r 2 - ri| = \J (x 2 - Xi) 2 + (y 2 - yi) 2 + (z 2 - z x ) 2 . 



It will be assumed that m\ is constant and that m 2 varies with respect the 
time. Taking into consideration this mass variation, Eqs. (la) and (lb) are 
written as 

d 2 v 1 _ Gm im2 

and 

d 2 r 2 Gm x m 2 . dr 2 /0 , 

where it has been defined m 2 as fn 2 = dm 2 /dt. Now, let us consider the usual 
relative, r, and center of mass, R, coordinates defined as 

r = r 2 — ri , and R = . (4) 

m 1 + m 2 

Let us then differentiate twice these coordinates with respect the time, taking 
into consideration the equations (2) and (3). Thus, the following equations 
are obtained 

r = - (mi +^g r -^r 2 (5) 

H m 2 

and 

- -m 2 . 2mim 2 . (m 1 + m 2 )m\ih 2 - 2m 1 rh 2 , 

R = r 2 + 7 -7T r H -r r . (6) 

mi + m 2 (mi + m 2 y (mi + m 2 y 

One sees that the relative motion does not decouple from the center of mass 
motion. So, these new coordinates are not really useful to deal with mass 
variation systems. In fact, using (4) , one has 

r 2 = R H ■ r , and r 2 = R H ■ r - ■ — r . (7) 

mi + m 2 mi + m 2 (mi + m 2 ) z 
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Substituting these expressions in (5) and (6), one can see more clearly this 
coupling, 



+ 



m\m\ 



m 2 (mi + m 2 ) 2 



m 2 
m 2 



R + 



mi 



mi + m 2 



and 



R = 



-m 2 



R + 



mim 2 



mi + m 2 (mi + m 2 ) 



(mi + m 2 )m 1 m 2 — m,im% 

-r H r . 

' (mi +m 2 ) 3 



(8) 



(9) 



However, one can consider the case for mi ^> m 2 (which is the case star-comet), 

— * 

and consider to put our reference system just on the first body (ri = 0). In 
this case, Eq. (3) becomes 



d 2 r 



Gmim 2 



m 2 r , 



where r = r 2 = (x,y,z). Using spherical coordinates (r,9,p), 

x = r sin 9 cos ip , y = r sin 9 sin ip , z = r cos 9 , 
Eq. (10) can be written as 



m 2 



dt 2 



Gmim 2 



+ m 2 f 



— iii 2 \r9 9 + r(p sin 9 (p 



where f, 9 and (p are unitary directional vectors, 

r = (sin 9 cos <p, sin 9 sin ip, cos 9) , with f=^6 , + (/9sin^(^ 



and 

= (— sine/?, cosy?, 0) , with </3 = — (p(sin9 r + cos# 9) . 

Since one has that r = rf, it follows that 

r = (f - r9 2 - rp sin 2 9)r + (2f 9 + r9 + r(p sin 9 cos 9)9 
+ (2rp sin # + rip sin 6* + 2r0^ cos 9)p , 

and Eq. (12) is discomposed in the following three equations 



m 2 (f — r9 — rip sin 2 9) = 



Gm 1 m 2 



m 2 r , 



m 2 (2r9 + r9 + rip sin 9 cos 9) = m 2 r9 , 



and 



(10) 



;n) 



(12) 



(13a) 



9 = (cos 6* cos v?, cos ^ sin (/?, — sin6>) , with 9 = — 9 r + ip cos 9 (p (136) 



(13c) 



(14) 

(15a) 
(156) 

m 2 (2rip sin 9 + rip sin 9 + 2rp>9 cos 9) = m 2 ripsin9 . (15c) 

Thus, one can restrict oneself to consider the case m 2 r m 0. For this case, it 
follows that ip = 0, and the resulting equations are 
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m 2 [r — rU ) = m 2 r , (16a) 

and 

m 2 {2r6 + r6) = . (166) 

Let m be the mass of the second body when this one is very far away from 
the first body (when a comet is very far away from the sun, the mass of the 
comet remains constant). Since m 2 7^ on (16b), the expression inside the 
parenthesis must be zero. In addition, one can multiply this expression by 
m r to get the following constant of motion 

P 9 = m o r 2 e . (17) 

Using this constant of motion in (16a), one obtains the equation [10] 

d 2 r _ Gm 1 | P$ m 2 (dr\ 
dt 2 r 2 m 2 r 3 m 2 \dt J 

This equation represents a dissipative system for m 2 > and anti-dissipative 
system for m 2 < 0. Suppose now that m 2 is a function of the distance between 
the first and second body, m 2 = m 2 (r). Therefore, it follows that 



dm 2 dm 2 dr 
dt dr dt 

and Eq. (18) can be written as 



(19) 



d 2 r _ Gmi | P[ m' 2 / dr\ 2 
dt 2 r 2 m 2 r 3 m 2 \dt J ' 

where m' 2 = dm 2 /dr. This equation can be seen as the following autonomous 
dynamical system [12] 

— = v — = Gmi + P g _ !Hk v 2 (21) 
dt dt r 2 m 2 r 3 m 2 

A constant of motion for this system is a function K = K(r,v) such that the 
following partial differential equation is satisfied [13] 



v 



OK , f-Gm 1 Pi m' 2 2 \dK 



) + ^P + ^fs-— v 2 ^ = 0. (22) 
or \ r z m z r A m 2 J ov 



This equation can be solved by the characteristic method [14] from which the 
following characteristic curve results 

C(r, v) = m\[r)v 2 + 2Gm x \ -^2 -JL / , (23 ) 



and the general solution of (22) is given by 

K(r,v) = F(C(r,v)) , (24) 
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where F is an arbitrary function of the characteristic curve. One can have a 
constant of motion with units of energy by selecting F as F = C/2m . That 
is, the constant of motion is given by 

K ( r v \ _ m 2( r ) v 2 ] Gm i f m l( r ) dr P e f m 2 (r) dr ^5) 
2m Q m Q J r 2 m 3 J r 3 



2. Lagrangian and Hamiltonian 

Given the time independent constant of motion (25), the Lagrangian of the 
system (20) can be obtained using the following known expression [11] 

r / % f K(r, v) dv . . 

L{r,v)=vJ . (26) 

Thus, the Lagrangian is given by 

L( r v ) = m 2( r ) t ,2 _ Gm i f m li. r ) dr | Pp f m l( r ) dr 
2m Q m Q J r 2 m 3 J r 3 

The generalized linear momentum (p = dL/dv) is 

P=^v, (28) 



and the Hamiltonian is 

m p 2 Gmi [ m 2 (r)dr P 2 r m 2 {r) dr 
H[r > P) ~ 2m 2 (r) + m n J r 2 m 3 J r 3 ' ^ 



2m\{r) 

Note from (25) and (29) that the constant of motion and Hamiltonian can be 
written as 

K(r,v) = ^±v> + V eff (r) (30) 







and 

m Q p 2 



mr ' p) = 2^) +V '" (r) - (31) 



where V e ff is the effective potential energy defined as 

Gmi f ml(r) dr Pg r m|(r) dr 
m 

This potential energy has an extreme value at the point 

p2 
r a 



Gmi r m 2 (r) dr P 2 r m 2 (r) dr 
Veffir) = ^7 J 'mil ■ (32Q) 



Gmim 2 



(326) 



which depends on m Q but it does not depend on the model for m 2 (r). One can 
see that this extreme value is a minimum for m 2 (r*) ^ 0, since one has that 

((fVefA _ (Gmim fm m 2 2 (r*) 

dr 2 )- R? >U - 
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On the other hand, because of the expression (28), one could expect different 
behavior of a trajectory in the phase space (r,v) and the phase space (r,p). 
The trajectory r{9) is found using the relation dr/dt = (dr/d0)0, and the Eq. 
(17) in (30) to get 



/ 

Jo, 



'dt>=*r m ' (T)dr , (34a) 
^ J„r^K-V e „(r) 



where K and Pg are determinate by the initial conditions, K = K(r ,v ) and 
Pg = m o r 2 o . The time of half of cycle of oscillation, Ty 2 , is obtained directly 
from Eq. (30) as 

fr-2 m2 (r) dr 
V2rn~oJri - V eff (r) 

where r\ and r 2 are the two return points deduced as the solution of the 
following equation 

V eff (n)=K, i = l,2. (34c) 



3. Model of Variable Mass 



As a possible application of (25) and (29), consider that a comet looses material 
as a result of the interaction with star wind in the following way (for one cycle 
of oscillation) 



TO 2 (r) 



m 00 \J\ — e ar incoming (v < 0) 

_ ?77je a(ri ~ r ) + m/(l — e~ ar ) outgoing (v > 0) 



(35) 



where ?77 00 or nif (where nif — 2rrii — m 00 by symmetry) is the mass of the 
comet very far away from the star (in each case), mi is the mass of the comet 
at the closets approach to the star (at a distance r x ), rrii = m 00 \J\ — e~ ari , and 
a is a factor that can be adjusted from experimental data. Thus, the effective 
potential (32a) has the following form for the incoming case (m = m 00 ) 



V 



(in) 



eff 



r) = 



Gm\m 



■{1-e 



+ 



P 2 



2m nn r 2 



[1 - e- ar ) 



+ 



Gm\m 00 a + 



a 2 P 2 



2m ( , 



Ei(-ar) + 



a Pie 



2„— ar 



2m nn r 



(36a) 



where Ei(x) is the exponential-integral function [15]. 
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For the outgoing case, one has m — rrif and 



Gmirrif 

r 2m fr 2 

Gmi(mie ari -m^ 2 



-2ar 



m f 

P$(mie ari -m f f 


r 

- e- 2ar 


m) 

+2Gm 1 (mie ari - mj 
2P|(m i e ari -m f ) 


2r 2 
r p~ ar 


) - — 

r 


m) 


[ 2r 2 ^ 



- 2aEi(-2ar) 

-2ar 



ae 



+ 2a 2 Ei(-2ar) 



— aEi(-ar) 



ae 



2r 



(366) 



where Pq is defined now as Pe = mf 2 Q. The extreme point of the effective 
potential (32b) for the incoming and outgoing cases is given by 



Gmim 2 



out G mi m 2 f 



(37) 



Given the definition (35), the constant of motion, Lagrangian, generalized 
linear momentum, and Hamiltonian are given by 



L (t) (r,v) 
■,(0 



p w (r, v) 



m|(r) 



and 



H®(r,p) 



m p z 



(38) 

(39) 
(40) 

(41) 



2m2(r) 

where i = in for the incoming case, and i = out for the outgoing case. As 
an example of illustration of this model, let us use the following parameters 
to estimate the dependence of several physical quantities with respect the 
parameter a, 

m 00 = l0 6 Kg ; P e = 10 17 Kg m 2 j sec ; and K = -8 x 10 23 J . (42) 

Fig. 1 shows the curves of V e ff(r) for several values of a (incoming case). As 
one can see from this figure, the location of the minimum does not change, but 
the minimum value of V e ff tends to disappear as a goes to zero. Fig. 2 shows 
the velocity (v) and normalized linear momentum (p/m ) as a function of r 
for several values of a and for the incoming case. All the trajectories start at 
r 2 = 200 and finish at ri(a). One can see the difference of the trajectories in 
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(a) with respect to (b) due to position dependence of the momentum, relation 
(40). 



4. Conclusions 

The Lagrangian, Hamiltonian and a constant of motion of the gravitational 
attraction of two bodies when one of them has variable mass were given. One 
feature of these quantities was the appearance of an effective potential, which 
is reduced (when m 2 = 0) to the usual gravitational effective potential of 
two bodies with fixed masses. Other feature was the distance dependence 
of the generalized linear momentum, Eq. (28). A model for comet-mass- 
variation was given which depends on the parameter a. A study was made of 
the dependence with respect to a of V e ff, minimum and maximum distance 
between the two bodies, and the trajectories in the spaces (r,v) and (r,p). 
These trajectories are qualitatively different for the same initial conditions 
due to the dependence of the linear momentum on the distance variable (r). 
Of course, the problem of the interaction comet-star with this model of the 
variation of mass deserves more complete analysis. The intention here with this 
example was to show explicitly the form of the constant of motion, Lagrangian, 
and Hamiltonian and to point out the different trajectories behavior in the 
spaces (r,v) and (r,p) arising from the constant of motion and Hamiltonian. 
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Figure Captions 

Fig. 1 v}jj\r) with the values of the parameters given on (42), for a — 1 
a = 0.01 (2); and a = 0.005 (3). 

Fig. 2 (a): Trajectories in the plane (r,v); (b): Trajectories in the plane (r 
a = 1 (1), a = 0.01 (2), and a = 0.005 (3). 
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